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ABSTRACT 

A recently developed constitutive model and a finite element formulation for pre- 
dicting the thermomechanical response of SMA hybrid composite (SMAHC) structures 
is briefly described. Attention is focused on constrained recovery behavior in this study, 
but the constitutive formulation is also capable of modeling restrained or free recovery. 
Numerical results are shown for glass/epoxy panel specimens with embedded Nitinol 
actuators subjected to thermal and acoustic loads. Control of thermal buckling, random 
response, sonic fatigue, and transmission loss are demonstrated and compared to conven- 
tional approaches including addition of conventional composite layers and a constrained 
layer damping treatment. Embedded SMA actuators are shown to be significantly more 
effective in dynamic response abatement applications than the conventional approaches 
and are attractive for combination with other passive and/or active approaches. 

INTRODUCTION 

Excitation levels for many aerospace vehicle structures are high due to engine 
noise and turbulent boundary layer fluctuating pressures. The ever-increasing need for 
weight-efficient structures leads to significant sonic fatigue and interior noise issues. 
Conventional aerospace structures typically employ passive treatments such as structural 
stiffening, constrained layer damping, and acoustic absorption materials to reduce the 
structural response and interior acoustic levels. These treatments typically suffer from 
substantial weight penalty and are often limited to relatively low temperatures and/or 
high frequencies. Smart materials may enable new structural design paradigms for high 
performance structures. There are many smart materials concepts, but some involving 
shape memory alloys (SMAs) appear to have significant advantages for applications 
involving adaptive stiffening and structural shape control. 

Shape memory alloys are a class of smart materials that exhibit a martensitic phase 
transformation when cooled through the transformation temperature range M s — Mj 
from the austenitic state, see figure 1. The shape memory effect (SME) can be described 
in simple terms in the following manner. A SMA can be easily deformed in the low 
temperature martensitic condition and can be returned to its original configuration by 
heating through the reverse transformation temperature range A s — Aj. This type of 
SME is termed free recovery. Conversely, in a constrained recovery configuration the 
SMA element is prevented from recovering the initial strain and a large tensile stress 
(recovery stress) is induced. A situation in which the specimen performs work (deforms 
under load) is called restrained recovery. Extensive work has been done to characterize 
these materials, both qualitatively through theoretical models [ 1—4] and quantitatively for 
particular alloys [5,6] . Compilations of papers have also been published, most recently 
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by Otsuka and Wayman [7J. Birman [8] gave a review of work done in the areas of alloy 
characterization, constitutive modeling, and applications. 

Although there are several alloy systems that exhibit the SME, Nitinol has the most 
engineering significance. The recovery stress (constrained recovery) for Nitinol is a 
function of initial strain (prestrain) and temperature, see figure 2 [5]. The yield strength 
of the austenitic phase is about ten times that of the corresponding martensitic phase. 
The austenitic Young’s modulus is also three to four times that of the martensite, see 
figure 3 [5J. 

Shape memory alloys have been investigated for a variety of applications since their 
discovery [9-12]. The concept of embedding SMA actuators in a composite laminate, 
termed a shape memory alloy hybrid composite (SMAHC), was introduced by Rogers 
and Robertshaw [13]. Subsequent studies have considered SMA actuators, external to 
the structure, for active vibration control [14-16] and shape control [17] of cantilevered 
beams. Other studies have proposed SMA actuators for vibration control of space struc- 
tures [18] and presented analytical formulations to demonstrate the vibration/structural 
acoustic control of SMAHC panels [19,20]. Studies of thermal post-buckling [21] and 
random dynamic response [22] suppression of SMAHC panels were presented using a 
thermoelastic FE formulation. A different approach was offered by Ro and Baz [23-25], 
where thermal, static, and dynamic analyses were developed for the case of SMA ac- 
tuators passing through sleeves in composite plates. Birman [26] investigated optimal 
distributions of SMA actuators in sleeves to improve the buckling performance of sand- 
wich panels under mechanical loading. 

SMAHC structures incorporating prestrained actuators are of particular interest in this 
study because of the potential for structural control due to tensile recovery stresses that are 
generated (constrained recovery) when the structure is placed in its service environment. 
In this way, an inherently elevated temperature (e.g. many aerospace applications) can 
be used to control structural response without control electronics or auxiliary power. 
The purpose of the present work is to demonstrate the dynamic response performance of 
SMAHC panels in comparison to conventional response abatement approaches. 

FORMULATION 
Constitutive Modeling 

Constitutive theories for SMAHC structures can be developed from the previously 
mentioned shape memory models [1-4], as was done by Boyd and Lagoudas [27] and 
Birman [28]. However, these models are difficult to use in practice. An alternative 
approach is to employ a constitutive model which makes use of experimental measurement 
of fundamental engineering properties. A new model of the latter type was recently 
developed by Turner [29]. This model casts the uniaxial thermoelastic constitutive relation 
for a SMA actuator, along the axis of the actuator, in terms of an effective coefficient 
of thermal expansion (CTE): 

= E a [e-i - J ai a (r)dr] (1) 

T 0 

where E a is the Young’s modulus of the SMA, is the mechanical strain in the 
1-direction, and o la is the “effective” (nonlinear) CTE. Note that this expression is 
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valid for constrained, restrained, or free recovery applications; only the empirical method 
of obtaining the thermal strain changes. A measure of the nonlinear CTE «1 a(T) 
over the temperature range of concern would be appropriate for free or restrained 
recovery applications. For constrained recovery applications, however, one must resort 
to measurement of the recovery stress and modulus. 

It can be shown that, for constrained recovery behavior, the nonlinear thermal strain 
in Equation (1) can be modeled by linear thermal expansion below the austenite start 
temperature and can be related to the SMA recovery stress oy and modulus by the 
equation 


T T 



at temperatures above A s . Note that in this case, the nonlinear thermoelastic nature of 
the SMA is still captured, albeit in a different way, because measurements of recovery 
stress and modulus versus temperature are inherently cumulative (integrated). Also, note 
that an experimental approach could be devised to measure appropriate recovery stress 
oy quantities for restrained recovery applications as an alternative to direct measurement 
of a la {T). The uniaxial SMA constitutive relation for the transverse direction has a form 
analogous to Equation (1). However, a- 2 a{T) is not related to the recovery stress, but is 
still nonlinear due to the differing martensitic and austenitic properties. 


Governing Equations 

The governing equations for the thermomechanical response of a structure subjected 
to combined steady-state thermal and out-of-plane dynamic mechanical loads can be 
derived through the use of a variational principle. The resulting finite element system of 
equations can be written in the following form [29J 
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or in compact notation as 

+ ([A-] - [A at] + ^[Al] + ^[A2]){A} = {F(t)} + {P AT } 


(3) 


(4) 


where [M] and [A’] are the usual system mass and linear stiffness matrices; [A’at] is the 
geometric stiffness matrix due to the thermal in-plane force vector {.Vat): [-V 1] and [A 2] 
are the first- and second-order nonlinear stiffness matrices which depend linearly and 
quadratically upon displacement {A}, respectively; { / ( / )} is the mechanical excitation 
load vector, and {Pat} is the thermal force vector. The subscripts b and m denote 
bending and membrane components, respectively, and the subscript B indicates that the 
corresponding stiffness matrix is due to the laminate bending-membrane coupling stiffness 
matrix [ B ]. Note that the stiffness matrices [A’], [A at], [Al], [A 2] and the thermal force 
vector {.Pat} are all temperature dependent. 


Proceedings of the 7th International Conference on Recent Advances in Structural Dynamics, The 
Institute of Sound and Vibration Research, University of Southampton, England, 24-27 July 2000. 



Three types of analyses are required to study the response of structures to thermal 
and dynamic mechanical loads, governed by Equations (4): (1) thermal buckling analysis, 
(2) thermal post-buckling analysis, and (3) dynamic analysis. Details of the constitutive 
model, FE formulation, and solution procedures can be found in Turner [29] . A Rayeigh’s 
integral approach was used to calculate radiated acoustic quantities and transmission loss, 
see Turner et al. [30]. 


RESULTS AND DISCUSSION 

The results presented in this study are intended to demonstrate the performance of 
SMAHC structures in comparison to conventional treatments. The results will focus on 
a single SMA actuator prestrain level of 5% because changes in SMA volume fraction 
have a more significant effect, the SMA volume fraction is more easily varied/controlled 
in fabrication, and the fatigue performance of Nitinol actuators degrades substantially at 
prestrain levels greater than 5%. The material system considered in this work consists 
of a E-glass/epoxy matrix material with Nitinol actuators embedded within the laminae 
along the direction of the glass fibers. The following material properties were used to 
generate the results presented in this section. 

Glass/Epoxy Nitinol 


Ei 

53.78 GPa 

(7,. 

from figure 2 

E, 

17.93 GPa 

E 

from figure 3 

G 12 

8.62 GPa 

V 

0.3 

”12 

0.25 

P 

6450 kg/rn 3 

P 

2031.8 kg/rn 3 

Ola 

6.61x10"° /°G, T<A S 

a i 

5.4x10-° /°G 


from equation 2, T>A S 

a 2 

30.6x10-° /°G 

a 2a 

6.61x10“° /°G, T<A S 


11. Ox 10" 6 /°C, T>A f 
interpolation, A s <T<Af 


Recall figure 2 characterizing the recovery stress for the Nitinol considered in this 
study [5]. The reverse transformation (austenitic) characteristic temperatures, under free 
recovery conditions, for the material represented in this plot are A s as 35°C (95°F) and 
Af «50°C (120°F). Note from figure 2 that Aj is increased by mechanical constraint 
and increasing prestrain, while A s is not significantly affected by either. Furthermore, 
the “knee” in the recovery stress curve becomes less distinct with increasing prestrain. 
These are well known phenomena and are attributable to stress increasingly inhibiting 
the formation of austenite and completion of the reverse transformation. These insights 
will be useful in later discussions. 

SMAHC Performance 

Consider a glass/epoxy panel with dimensions of 0.36x0.25x0.001 m (14x10x0.045 
in.) and (45/0/-45/90/0/90/-45/0/45) lamination. The reason for the middle 0°-layer will 
be made clear subsequently. Static and dynamic performance results for a conventional 
(baseline) and two SMAHC laminates with clamped boundary conditions are shown in 
table 1 for ambient 21°C (70°F) and elevated temperature 82°C (180°F) conditions. The 
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SMAHC laminates in this table have a SMA volume fraction of 22%. The dynamic 
response is forced by normally-incident, random acoustic pressure with an overall sound 
pressure level (OASPL) of 120 dB (ref 20 pPa) and a bandwidth of 1 0-400 Hz. A 
uniform modal critical damping ratio of 1% was used in the dynamic calculations. 

Some general comments about this data are as follows. The buckling and resonance 
behavior of the SMAHC panels relative to baseline panel at ambient temperature is at- 
tributable to differences in material properties, primarily in the 2-direction in principal 
material coordinates. Discussion of these differences is omitted here for brevity. Com- 
parison of the two SMAHC panels, one with SMA in all layers and the other containing 
SMA in the outer 0° layers and in the 90° layers, indicates some small performance gain 
is achieved by concentrating the embedded SMA in only a few layers. This approach is 
attractive from a fabrication point of view, particularly for automated fabrication where it 
is easy to imagine replacing strategic spools of prepreg with SMA ribbon in a continuous 
consolidation process. The middle 0° -layer is necessary to separate the two 90° layers 
containing a substantial SMA volume fraction, thereby avoiding consolidation problems 
during cure. The remaining results/discussions will focus on SMAHC laminates with 
embedded SMA concentrated in the outer 0° layers and in the 90° layers. 

It can be seen that activation of the embedded SMA has strong effects on the 
critical buckling temperatures, resonance frequencies, and RMS dynamic responses. The 
maximum displacement response power spectral density (PSD) of the SMAHC panel is 
shown in comparison to that of the baseline panel in figure 4. It is clear that the embedded 
SMA has little effect at ambient temperature. At 82 (180°F, AT=61°C), however, 
substantial benefit from the SMA recovery stress is evident in shifting the fundamental 
frequency, shifting other modes out of bandwidth, and reducing peak responses. Thus, 
performance can be enhanced by shifting resonance out of the excitation bandwidth and 
by response reduction within the excitation bandwidth. 

The data in table 1 also indicates that SMAHC panels can exhibit two critical 
buckling temperatures (A T cr for inactive and active SMA states), as was previously 
reported [22,29]. This occurs for cases in which A s is greater than ambient temperature. 
This phenomenon is exemplified in figure 5, where the normalized, maximum post- 
buckling deflection versus temperature increase is shown for SMAHC laminates with 
varying SMA volume fraction. It can be seen that a volume fraction of 4% is insufficient 
to eliminate thermal post-buckling, but its effect in reducing the deflection is clearly 
observed. As little as 13% SMA will eliminate thermal post-buckling over a small 
range of temperatures; higher concentrations of SMA (e.g. 22%) can greatly extend 
the temperature range. This behavior can be explained as follows. Thermal expansion 
may induce buckling at temperatures below ,l, s . The SMA is activated at A s and the 
recovery stress rapidly increases with temperature, which overcomes the compressive 
thermal stress and reduces the thermal deflection (potentially rendering the structure flat 
again). At higher temperatures, the rate of increase of the recovery stress with temperature 
diminishes and may be surpassed by the rate of increase in matrix compressive thermal 
stress. Further increases in temperature will cause net-compressive stress development, 
eventually leading to a second incidence of thermal buckling. Note that this effect can 
be easily avoided in practice by selecting an alloy with slightly higher Nickel content 
[9], which can lower A s to ambient temperature or below and result in a single, high- 
temperature buckling event. 
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The thermal post-buckling observations have strong implications for the panels’ 
dynamic response. An example is shown in figure 6, where the maximum displacement 
PSD of the SMAHC panel (r, ; =22%) is shown at ambient and two elevated temperatures. 
At 82°C (180°F, AT=61°C), substantial benefit of the embedded SMA is evident, while 
further heating to 121°C (250°F, AT=100°C) greatly diminishes the performance (RMS 
W max =0.0927e-3 m) because continued thermal expansion has partially overcome the 
beneficial effects of the recovery stress. This particular material system, glass/epoxy 
matrix with 22% Nitinol volume fraction at 5% prestrain, achieves its “optimal” dynamic 
response at approximately 82° C (180°F), which corresponds to the “knee” in the recovery 
stress curve, see figure 2. This result also holds for analogous SMAHC panels (V,,=22%) 
with embedded Nitinol having 3% and 4% prestrain, which exhibit “optimal” dynamic 
performance at their respective “knee” temperatures of 54°C (130°F) and 66°C (150°F). 
The optimal temperature for dynamic performance will be designated by A 0 and is 
governed by the same thermoelastic balance as that described in the thermal post-buckling 
discussion. In general, with a different material system or at higher prestrain levels (in the 
absence of a distinct “knee”), A 0 must be determined numerically or by comparing the 
slope of the recovery stress curve with that of the matrix thermal stress versus temperature. 

Selection of an alloy composition and processing combination that will place A 0 near 
the application temperature is obviously recommended for dynamic control applications. 
An alternate approach to achieve similar performance would entail increasing the SMA 
volume fraction and/or prestrain, which will result in a less efficient and/or durable design. 
The above discussions imply potential for a delicate balance between avoiding thermal 
buckling by imposing A s to be less than ambient temperature and “optimizing” dynamic 
performance by specifying A 0 to be near the application temperature. Thermal buckling 
can also be avoided by boundary expansion, a common approach for aerospace structures, 
so A s can be neglected in such cases and A 0 should be used to guide efficient design. 

SMAHC vs. Conventional 

The performance of the SMAHC laminate will now be compared with conventional 
approaches including increased stiffness through additional composite layers and a con- 
strained layer damping (CLD) treatment. The conventional approaches were analyzed at 
ambient temperature to simulate the conditions of boundary expansion, i.e. no thermal 
buckling. No attempt was made to account for changes in matrix properties with temper- 
ature. A variety of dynamic response performance characteristics are compiled in table 2 
for convenience of comparison. The SMAHC panel achieves a 77% RMS displacement 
response reduction, relative to the baseline 9-layer glass/epoxy panel, with a 48% mass 
increase. A conventional composite laminate of 13 layers (45/0/-45/90/45/0/90/0/45/90/- 
45/0/45) exhibits a similar mass increase (44%), but only achieves a 61% RMS displace- 
ment reduction. A RMS displacement response of 0.0474e-3 m (nearly equivalent to that 
of the SMAHC panel) can be achieved with a 16-layer (45/0/-45/90)2 S conventional lam- 
inate, but at the expense of a 78% mass increase. A typical CLD treatment was modeled 
as a 19.16 N/m 2 (0.2 lbf/ft 2 ) weight addition and a uniform modal critical damping ratio 
of 6%, which are considered to be representative values. This treatment also happens to 
result in a mass increase (42%) similar to that of the SMAHC, but the RMS displacement 
response reduction of 63% is again inferior to the SMAHC. Additional insights are more 
clear from comparisons of the displacement response spectra. 
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A comparison of the respective maximum displacement response PSDs is shown in 
figure 7. As anticipated, the means by which the three systems achieve dynamic response 
reduction is quite different. The SMAHC achieves a much more weight-efficient stiffen- 
ing effect relative to that of adding layers, while the CLD treatment achieves reduction 
primarily by damping enhancement. The weight-efficient stiffening effect of the em- 
bedded SMA is advantageous for combination with other passive treatments, which are 
typically more effective at higher frequencies. SMAHCs are also attractive for combi- 
nation with control approaches involving piezoelectric actuators, which typically exhibit 
actuation authority roll-off below 200 Hz. Such combined active-passive approaches 
could achieve enormous dynamic response reductions by using the embedded SMA to 
shift resonance frequencies and provide peak response control by stiffness enhancement, 
while the active approach would provide enhanced damping at resonances still within 
the excitation bandwidth. 

The dynamic response results have direct implication for sonic fatigue life prediction. 
Consider an analogous acoustic pressure loading condition with an OASPL of 150 dB (ref 
20/ /Pa). It is clear that linear dynamic response analysis is no longer valid, but its use at 
high response levels is still common practice in the aerospace industry for sonic fatigue 
design. This condition is more representative of a sonic fatigue environment and will 
produce more realistic life predictions, while it is understood that conservative estimates 
will be obtained. A representative strain-life (S-N) relationship for glass/epoxy laminates 
was modeled by e = 7385e -6 (iV) - ' and the maximum RMS normal strain was 
taken as the failure criterion. In all cases, the maximum strain was the transverse normal 
strain at the mid-span of the long boundary. Under these conditions, the baseline panel 
has an expected life of 63.7k cycles. The corresponding life predictions for the SMAHC, 
13-layer conventional, and CLD panels are 205M, 2.84M, and 37. 8M cycles, respectively. 
These estimates are in direct relation with the corresponding peak displacements in figure 
7. It is intuitive that measures to further attenuate the SMAHC peak response would result 
in excellent sonic fatigue performance. The previously mentioned active-passive system, 
for example, might be a particularly attractive solution. 

A plot of the transmission loss (TL) versus frequency for the SMAHC panel at two 
temperatures is shown in comparison to the baseline panel in figure 8. The main effect 
of the embedded SMA in its most effective state at 82°C (180°F) is a vast increase 
in the stiffness-controlled portion of the TL. This improvement is diminished at higher 
temperatures (e.g. 121°C) for the same reasons as stated above for the dynamic response. 
Although the stiffness-controlled TL for the SMAHC panel shows dramatic improvement, 
the TL over the 400 Hz bandwidth is actually diminished because the mass-law effect 
has been moved to higher frequencies. The TL for the SMAHC panel at 82°C (180°F) 
is shown in comparison to that for the two conventional approaches in figure 9. It can 
be seen that the SMAHC panel is superior over a 200 Hz bandwidth, but not as effective 
at mid-range (resonance-controlled) frequencies due to the shift in mass-law effects. 
Combination of the SMAHC with other passive and/or active approaches is again an 
attractive alternative for this frequency range. At high frequencies (mass controlled) the 
three approaches will produce nearly equivalent TL. 

CONCLUSIONS 

A recently developed thermoelastic constitutive model and a finite element formula- 
tion for predicting the thermomechanical response of shape memory alloy (SMA) hybrid 
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composite structures (SMAHC) have been briefly described. Constrained recovery behav- 
ior of embedded SMA actuators is the focus of this study, but the constitutive formulation 
is also capable of modeling restrained or free recovery behavior. The constitutive model 
captures the material nonlinearity of the SMA with temperature. The constitutive and 
finite element formulations are in a form that is amenable to implementation in a com- 
mercial finite element code. 

Results are shown for SMAHC panels consisting of a glass/epoxy matrix material 
system with embedded Nitinol actuators having a prestrain of 5% subjected to uniform 
thermal and acoustic pressure loads. Performance gains can be achieved by concentrating 
the SMA actuators in specific layers, which is also driven by fabrication considerations. 
SMAHC panels can exhibit two critical buckling temperatures for A s > T 0 , which can be 
avoided by proper alloy composition selection. SMAHC structures exhibit an “optimal” 
operating temperature (A,,) for dynamic response performance. A 0 is a function of 
the matrix material and SMA composition, volume fraction, and prestrain level. The 
implied conflict between specifying A s < T 0 to avoid thermal buckling and A 0 ps T 
to “optimize” dynamic response can be resolved by taking the conventional approach of 
allowing boundary expansion to prevent buckling and using A 0 % T to guide structural 
design. 

The performance of a SMAHC panel was compared to that of a baseline panel 
and two conventional approaches including additional conventional composite layers 
and a constrained layer damping (CLD) treatment. The three response abatement 
approaches have similar added mass relative to the baseline panel, but the SMAHC 
panel exhibits superior dynamic response and sonic fatigue behavior. The SMAHC panel 
also exhibits superior transmission loss (TL) characteristics in the low frequency, stiffness 
controlled region; up to 200 Hz in this case. The SMAHC panel is less effective for 
TL at higher frequencies due to the shift in mass-law effects. However, the weight- 
efficient stiffening effect of the embedded SMA makes it ideal for low frequency noise 
transmission abatement and for combination with other passive or active approaches, 
which are typically more effective at higher frequencies. These combined approaches 
have significant implications for dynamic response/sonic fatigue also. 
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Table 1: Results for conventional and SMAHC panels. 


SMA Layer 
Content 

T=21°C, Inactive SMA 

T=82°C, Active SMA 

AT cr , °C 

ft, Hz 

RMS 

W max? m 

AT cr , °C 

ft, Hz 

RMS 

Wmaxj m 

None 

4.99 

87.48 

0.2062e-.3 

4.99 

— 

— 

All 

5.29 

71.11 

0. 1899e-3 

113.6 

173.5 

0.0482e-3 

Outer 0°, 90° 

5.43 

87.57 

0.2058e-3 

122.1 

226.4 

0.0476e-3 


Table 2: Results from baseline, SMAHC, and panels 
with conventional response abatement treatments. 


Case 

Added 

Mass 

ft, 

Hz 

RMS 

W max , m 

TL, 

dB 

Max RMS 
pe ® 120dB 

Life ® 
150 dB 

Baseline 

0% 

87.48 

0.2062e-.3 

17.9 

42.0 

63.7e+3 

SMAHC ® 
82°C 

48% 

226.4 

0.0476e-3 

14.3 

12.0 

205e+6 

13- Layer 

44% 

128.3 

0.0805e-3 

19.1 

23.3 

2.84e+6 

16- Layer 

78% 

158.6 

0.0474e-3 

20.2 

16.8 

2.3.4e+6 

Baseline-CLD 

42% 

73.4 

0.0768e-3 

27.4 

15.6 

37.8e+6 
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Figure 2: Nitinol recovery stress versus temperature and initial strain [5]. 



Temperature, °C 


Figure 3: Nitinol Young’s modulus versus temperature [5]. 


Proceedings of the 7th International Conference on Recent Advances in Structural Dynamics, The 
Institute of Sound and Vibration Research, University of Southampton, England, 24-27 July 2000. 







Figure 4: W max PSD for the baseline and SMAHC laminates. 



Figure 5: Post-buckling deflection of SMAHC laminates versus AT and v a . 



Figure 6: W max PSD for a SMAHC laminate at three temperatures. 
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Figure 7: II .. . r : l’SI) comparison of two SMAHCs with conventional treatments. 



Figure 8: TL versus frequency for the baseline and SMAHC laminates. 



Figure 9: TL comparison of two SMAHCs with conventional treatments. 
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